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Abstract 

Background: MicroRNAs (miRNAs) are small, endogenously expressed, non-coding RNA molecules involved in 
gene transcription and expression that combine with specific mRNA site of target genes to inhibit protein synthesis 
or degrade mRNA. Since the first plant miRNA was reported in 2002, numerous new miRNAs and their targets have 
been discovered via high-throughput sequencing and computational approaches. However, the genetic variation of 
miRNA genes is poorly understood due to the lack of miRNA-specific DNA markers. 

Results: To study the genetic variation and map miRNAs and their putative target genes in cotton, we designed 
specific primers based on pre-miRNAs and published putative target genes. A total of 83 pre-miRNA primers and 
1,255 putative target gene primers were surveyed, and 9 pre-miRNA polymorphic loci were mapped on 7 of the 26 
tetraploid cotton chromosomes. Furthermore, 156 polymorphic loci of the target genes were mapped on the 
cotton genome. To map more miRNA loci, miRNA-based SRAP (sequence-related amplified polymorphism) markers 
were used to map an additional 54 polymorphic loci on the cotton genome with the exception of ChrOl, Chr22, 
and Chr24. Finally, a network between miRNAs and their targets was constructed. All pre-miRNAs and 98 putative 
target genes were selected for RT-PCR analysis, revealing unique expression patterns across different fiber 
development stages between the mapping parents. 

Conclusions: Our data provide an overview of miRNAs, their putative targets, and their network in cotton as well 
as comparative expression analyses between Gossypium hirsutum and G. barbadense. These data provide a 
foundation for understanding miRNA regulation during cotton fiber development. 
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Background 

MicroRNAs (miRNAs) are a class of non-coding small 
RNA molecules, generally 21 nucleotides (nt) in length 
that regulate critical functions in plant and animal de- 
velopment. miRNAs arise from an -yO-nt stem-loop 
RNA precursor that is cleaved and modified by Ribo- 
nuclease III to create an ~21-nt single-stranded miRNA 
[1,2]. Generally, every mature miRNA derives from one 
of the precursor s arms and all are located in intergenic 
regions, indicating that miRNA transcription is inde- 
pendent of other genes and has its own transcriptional 
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regulatory mechanisms. miRNAs regulate gene expres- 
sion at the post-transcriptional level through mRNA 
cleavage or translational inhibition in both animals and 
plants [3-7]. In plant cells, the miRNA sequence is 
almost completely complementary with target gene 
mRNA, through a mechanism similar to RNA interfer- 
ence (RNAi), leading to the degradation of target gene 
mRNA [8]. Recently, numerous miRNAs have been 
reported through cloning, high-throughput small RNA 
sequencing, and computational approaches based on se- 
quence similarities and secondary structure predictions 
[4]. Because plant miRNAs are almost exactly comple- 
mentary to their corresponding target sequences [5], 
target gene prediction is relatively straightforward and 
mainly through bioinformatics method. 
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In 1993, the first miRNA {lin-4 gene) was reported in 
Caenorhabditis elegance [9]. Then, Reinhart et al [10] 
found the second metachronous switch gene, let-7, in C. 
elegans. Recently, with bioinformatics, molecular cloning 
technology improvements, and the establishment of model 
species cDNA libraries, hundreds of small RNA molecules 
have been successively reported in C. elegans [11,12], 
Arabidopsis [13,14], rice [15], maize [16], wheat [17], moss 
[18] and in cotton [19-22]. Almost 3,600 miRNAs have 
been annotated in Magnoliophyta including 291 from 
Arabidopsis thaliana and 40 from Gossypium species, as of 
April 21, 2012 (http://www.mirbase.org/cgi-bin/browse.pl). 
Recently, with the sequencing of the G. raimondii (D5) gen- 
ome, many miRNA precursors and miRNAs were identified 
[23-25]. However, their distribution on the tetraploid cotton 
genome is unknown. To date, few miRNAs have been func- 
tionally characterized in cotton, and almost no predicted 
target genes have been verified by experiments in cotton. 

Cotton is not only a fundamental world commodity 
that provides an important natural material, it is also 
an important oil crop. Since the first miRNA was re- 
ported, studies suggest that plant miRNAs negatively 
regulate target genes involved in plant development, 
organ morphogenesis, auxin signaling, and environ- 
mental stress responses [14,26-30]. MiRNAs may also 
play an important role in cotton fiber development and 
the study of miRNAs has far-reaching significance to 
improve fiber quality and yield [31]. 

At this time, most miRNAs and target gene research is 
focused on functional verification through gene over- 
expression, gene interference, and related methods, but 
there is paucity in development of miRNA markers and 
their relationship with phenotypes. Only one report in 
the literature described miRNA markers- miRNA- AFLP 
(amplified fragment length polymorphism) [32], but no 
miRNAs or their target genes have been genetically 
mapped in cotton. 

To address this gap and to map miRNAs to under- 
stand their genomic distribution in cotton, we down- 
loaded pre-miRNA sequences reported in the literature 
and within the miRBase database (Additional file 1: 
Table SI). Putative target genes were predicted by 
psRNATarget (http://plantgrn.noble.org/psRNATarget). Spe- 
cific primers were then designed based on the pre-miRNA 
sequences and the target sites of the target genes, and 
miRNA-based SRAP marker analysis was performed. 
Primers were also designed surrounding the target sites in 
the putative target genes to indicate their chromosomal 
distribution. Then, a network diagram between miRNAs 
and their putative targets was constructed to reveal 
regulatory relationships. Finally, RT-PCR analysis was 
conducted to detect expression differences between G. 
hirsutum and G barbadense in selected pre-miRNAs 
and putative target genes. 



Results and discussion 

Strategies of development markers for miRNAs and their 
target genes 

Primers were designed for 123 pre-miRNA sequences and 
identical primers were eliminated with a nucleotide BLAST 
search. Finally, 83 pairs of pre-miRNA primers were 
obtained (Additional file 2: Table S2). Then, these primers 
were screened using single-stranded conformation poly- 
morphism (SSCP), of which 12 primers yielding 13 poly- 
morphic loci. The rate of primer polymorphism was 14.5%. 

A total of 1,255 primer pairs were obtained for the 
1,399 target genes of the miRNAs after eliminating 
identical primers as described above (Additional file 3: 
Table S3). Of these, 147 primer pairs produced 161 
polymorphic loci after SSCP analysis. The rate of pri- 
mer polymorphism was 11.6%. 

Based on miRNA- AFLP technology, the miRNA-SRAP 
marker technique was performed to enrich miRNA 
markers. Among the 2,944 miRNA-SRAP primer combi- 
nations, ~8 scorable bands per primer combination were 
typically produced. Exceptions were miR319-SRAP pri- 
mer combinations which generated only 0-2 scorable 
bands. These results are consistent with those obtained 
from miRNA-based AFLP marker techniques [32]. Fifty- 
five primers produced 59 polymorphic loci (Additional 
file 4: Table S4), and the rate of primer polymorphism 
was 1.87%. No polymorphic amplification products were 
obtained for miR156, miR319, miR390, miR398, and 
miR399 when combined with SRAP primers. 

Pang et al. [32] explored the possibility of using miRNA 
genes as markers and identified DNA sequences of poten- 
tial pre-miRNA. miRNA- AFLP analysis provides a reliable 
targeted genotyping strategy to assess genetic diversity 
among cotton species. The SRAP primers can target exons 
and introns respectively and generate polymorphism based 
them [33]; they also can combine with primers desighed 
from candidate gene regions, which is also called TRAP 
(Target Region Amplification Polymorphism) [34]. There- 
fore, when primers designed from conservative miRNAs 
and their complementary sequences combining with SRAP 
primers, miRNAs or their flanking sequences can be amp- 
lified. This is a highly effective strategy to examine miRNA 
gene distribution as well as more feasible, simple, and effi- 
cient than the miRNA-AFLP technique with respect to 
typing miRNA markers. Also, miRNA-SRAP markers may 
be functional genes; therefore, the miRNA-SRAP tech- 
nique has important significance to the study of miRNAs. 

Verification of miRNA-SRAP marker techniques 

To verify whether the miRNA-SRAP marker technique 
reliably amplified miRNA sequences, PGR products were 
cloned and sequenced from a sample of 7 random poly- 
morphic miRNA-SRAP primer combinations. Then, 
miRNA primers and SRAP primers were identified as 
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depicted in Figure 1. The figure depicts reliably ampli- 
fied products using the primer combinations instead of 
amplified products from individual miRNA primers or 
SRAP primers. 

Even distribution of miRNA markers and preferential 
distribution of target gene markers in the cotton genome 

After linkage analysis, 9 of the 13 polymorphic loci of pre- 
miRNAs were mapped on 7 cotton chromosomes (ChrOl, 
Chr04, Chris, ChrlS, Chrl6, Chr23 and Chr26) based on 
SSCP analysis which can discover minor sequence muta- 
tions and reveal more polymorphisms to map genes [35]. 
For miRNA-based SRAP markers, 54 of the 59 poly- 
morphic loci were mapped on 23 cotton chromosomes, ex- 
cept for ChrOl, Chr22, and Chr24. These data suggest that 
miRNA markers are widely distributed on 24 chromo- 
somes, except for Chr22 and Chr24, for which no miRNA 
markers were identified. Apparently, relatively more loci oc- 
curred on ChrlO (5 loci) and Chrl9 (6 loci). Thirty- two loci 
were mapped on 13 chromosomes of the At sub-genome, 
and 31 loci were mapped on 11 chromosomes of the Dt 
sub-genome (Additional file 5: Figure SI, Additional file 6: 
Table S5). In conclusion, these 63 miRNA loci were evenly 
distributed on 24 chromosomes, except Chr22 and Chr24. 

Of the 161 polymorphic loci for the putative miRNA 
target genes, 156 were mapped on the cotton genome: 
59 loci were mapped on 13 chromosomes of the At 
sub-genome, and 97 loci were mapped on 13 chromo- 
somes of the Dt sub-genome (Additional file 5: Figure SI, 
Additional file 6: Table S5). Overall, unlike miRNA 
markers, target gene markers were preferentially distributed 



on the Dt sub-genome. Specifically, chromosomes Chrl3, 
Chris, Chrl9 and Chr21 contained a total of 47 target 
gene marker loci, i.e., 30% of the total target gene marker 
loci. Comparatively fewer less loci were on ChrOl, Chr03, 
Chr06, Chr07, ChrOS and Chr09. Surprisingly, 95 of the 
156 target gene polymorphic loci belonged to target genes 
from the miR414 family, whereas 61 loci belonged to 26 
other miRNAs families. 

To date, no miRNA genes and their target genes have 
been genetically mapped in cotton. Pang et al. [32] de- 
veloped conserved miRNA gene markers and explored 
the possibility of an miRNA-AFLP marker technique. 
Then, the amplified products were sequenced and com- 
pared for homology using cotton expressed sequence 
tags. However, genetic mapping of miRNA markers has 
not been reported. Thus, miRNAs and their target gene 
markers reported here were genetically mapped in our 
interspecific G. hirsutum x G. barbadense population. 
Such genetic mapping of a segregating population pro- 
vided an overview of the chromosomal distribution of 
miRNAs and their target genes in cotton. 

Network diagram revealed by relationships between 
miRNAs and their target genes 

In order to excavate the valuable miRNAs and putative 
targets for further research, miRNAs and putative target 
genes, mapped on the interspecific BCi genetic linkage 
map, were used to create a network diagram based on 
their regulatory relationships (Figure 2, Table 1) by Circos 
(http://circos.ca/). Table 1 depicts 39 miRNAs and 24 tar- 
get genes that are contained in this network diagram. 



miR165'Me21 



miR16S'Me22 



miR164-Me40 



miR16E-IVlel7 



miR403-Em7 



miR164'Em33 



miR169a'Me27 



TGAGTCCAAACCGGAGTGGAGTTGATCGGGGAAAAAATTAAAATTATGGATTAGCCCT 
TGAAGAACTCGCmACnATTTTGAACAGTTAAGTTCATATAGAACmCCACTTTATTT 
TATTTTATGT AGGTTTGAAGCCreGrCCGA 

TGAGTCCAAACCGGAGGAGACCAAGGAGTTGCAAAGACAAGTCAATGAACTGATG6A 
AAAAGGCTACATTCGCGAGAGCTTAAGTCCCTGTGCTATTCCCGAACGACTCTCTCGAAC 
GACGATCGGGCAACTAAGTGAATTTGACTCATCA rCCCGArcraCACC/WGCG 

TGAGTCCAAACCGGACCATACCATAGTTCCTCCAACCTATCGATTCAGTTTCATATCAATT 
CATGCATTCACATTCCATATACTCATCAATAGGACATTGAAGCATTTTCATAAATCAATTT 
ATmCGnACAATCAAACGACTAAACGGCCmCACCCATCCTACGATAAATTTTAfGOl 
CSmCCTGCTTCrCCA 

TGAGTCCAAACCGGATTGAACTAGAAnGCCATCATCCGAACrrGCCAAAGATTGAAAT 
TTGTCTCACCATCGAACTTCTCAATnCAAACAnGTTATTTCCATATCTGAATAGACTGAT 
CTTTGAACTAGCTCTGATACCACTTGTTAGGATTG7CCCGA TC7GCACCAAGCGjil 

GACTGCGTACGAATTATGGCGAATGACGTAAAAAGTAACGCACCGGCTCCCGCTGAAG 
GGGCAATGCCATCTGATAATCAACCTCCCACAGTTAGCCAGGGAAGAGGCAGrGCAr 

GGAGAAGCAGGGCACGraCACATCGAAGGTGTTAAGAATGGCATGGAAGTGAACGTT 
CAAAGGGCGAGGTCGATGGATGTAGAACTACCCAAAACTTGAAAACTTTCTGAGTGATA 
GAGTACATCAGCGGTCAGTCAAAAATTCGTACQCAGTC 

AGCCAAGGArG^CTTGCCGAGAAGGTGAAGGGAAAGTGAAATAAGGGAACCAATGAA 
AATCGCTGAAAACTGACAACATGCATACCGGTTTGGAaCA 



Figure 1 Verification of miRNA-based SRAP marker technology by sequencing. SRAP primers are in bold. miRNA primers are in Italic, 
bolded, and in shadow. 
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Figure 2 Network diagram revealed by relationships between miRNA families and their target genes. Red circle: MiRNAs; black circle: 
Target genes. The scale marked on each chromosome represents genetic map distance (cM). miRNA primers and target primers are marked on 
the chromosome, and the miRNA family and target gene relationships are indicted by each color. 



Some miRNAs were associated with more than one target 
gene, and some target genes were regulated by more than 
one miRNA family. Individual miRNAs and their target 
genes are clearly and intuitively depicted in the network 
diagram (Figure 2, Additional file 7: Figure S2). As shown 
in Additional file 7: Figure S2, miR393 and its target 
TC38219 (Additional file 7: Figure S2I), and miR172 and 
its target TC30742 were located on the same chromosome 
(Additional file 7: Figure S2H). Also, miR393 and its target 



TC38219 were proximal on the same chromosome, sepa- 
rated only by 1.78 cM; whereas, miR172-MelO and its tar- 
get TC30742 were separated by 30.72 cM. miR172-Me49b 
and its target TC30742 were separated by 11.41 cM 
(Table 1). Additional file 7H shows that miR172 and its 
target TC30742 (both located on ChrlO), and another 
target TC41731 (located on Chr20), are found on 
homologous chromosomes (ChrlO and Chr20). The 
miR396 (Additional file 7: Figure S2K) and miR396a 
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Table 1 miRNAs and their target genes localization on the interspecific BCi 


linkage map 






miRNAs name 


Chromosome 


Subgenome 


Genetic map 


Target name 


Chromosome 


Subgenome 


Genetic map 








distance (cM) 








distance (cM) 


miR156c-DX535560* 


Chr23 


D 


154.258 


T-TC29212 


Chrll 


A 


112.712 










T-CO079328 


Chr21 


D 


56.729 










T-TC37339 


Chr21 


D 


220.676 


miR156c-DX535560* 


Chr23 


D 


154.258 


T-TC35884A 


Chrll 


A 


214.024 










T-TC36356A 


Chrl5 


D 


1 18.218 


miR157b-CO076888* 


Chr23 


D 


150.354 


T-TC35884A 


Chrl 1 


A 


214.024 










T-TC36356A 


Chrl5 


D 


118.218 


miR159-Me29 


Chr05 


A 


165.2 


T-TC32762 


Chrll 


A 


94.707 


miR159-Me24 


Chrl2 


A 


156.471 


T-TC32762 


Chrl 1 


A 


94.707 


nniR164-Me4a* 


Chrll 


A 


146.281 


T-TC41649 


Chr05 


A 


70.705 


miR164-Me4b* 


Chr20 


D 


26.305 


T-TC41649 


Chr05 


A 


70.705 










T-CO 108738 


Chr09 


A 


124.81 1 


miR164-Mel2* 


Chrl4 


D 


157.964 


T-TC41649 


Chr05 


A 


70.705 










T-CO 108738 


Chr09 


A 


124.811 


miR164-Me29* 


Chrl7 


D 


79.157 


T-TC41649 


Chr05 


A 


70.705 










T-CO 108738 


Chr09 


A 


124.811 


miR164-Me40* 


Chrl3 


A 


91.878 


T-TC41649 


Chr05 


A 


70.705 










T-CO 108738 


Chr09 


A 


124.81 1 


miR164-Em33* 


Chrl7 


D 


82.728 


T-TC41649 


Chr05 


A 


70.705 










T-CO 108738 


Chr09 


A 


124.81 1 


miR164-Em20* 


Chr26 


D 


143.941 


T-TC41649 


Chr05 


A 


70.705 










T-CO 108738 


Chr09 


A 


124.811 


miR167-Me7* 


Chr25 


D 


126.296 


T-TC31104 


Chrl2 


A 


169.342 










T-TC38767 


Chr23 


D 


25.125 


miR167-Me49a* 


Chr03 


A 


67.893 


T-TC31104 


Chrl2 


A 


169.342 










T-TC38767 


Chr23 


D 


25.125 


miR167-Me49b* 


Chr05 


A 


198.84 


T-TC31104 


Chrl2 


A 


169.342 










T-TC38767 


Chr23 


D 


25.125 


nniR169a-Me27 


Chrl9 


D 


18.195 


T-TC29763A 


Chrl3 


A 


153.525 


miR169a-Em5 


Chrl8 


D 


84.73 


T-TC29763A 


Chrl3 


A 


153.525 


nniR169a-EnnlO 


Chr03 


A 


75.625 


T-TC29763A 


Chrl3 


A 


153.525 


nniR169a-Eml2 


Chris 


D 


134.58 


T-TC29763A 


Chrl3 


A 


153.525 


miR169b-Me45 


ChrOS 


A 


131.477 


T-TC29763A 


Chrl3 


A 


153.525 


nniR169b-EnnlO 


Chrl9 


D 


64.043 


T-TC29763 A 


Chrl3 


A 


153.525 


miR169b-Em29b 


ChrlO 


A 


87.599 


T-TC29763A 


Lnr 1 i 


A 

A 


1 53.525 


miR172-Me6* 


Chr21 


D 


149.272 


T-TC30742 


ChrlO 


A 


123.541 










T-TC33909 


Chrl6 


D 


82.967 










T-TC41731 


Chr20 


D 


80.438 


miR172-MelO* 


ChrlO 


A 


92.823 


T-TC30742 


ChrlO 


A 


123.541 










T-TC33909 


Chrl6 


D 


82.967 










T-TC41731 


Chr20 


D 


80.438 
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Table 1 miRNAs and their target genes localization on the interspecific BCi linkage map (Continued) 



miR172-Me49a* 


Chrl3 


A 


76.047 


T-TC30742 


ChrlO 


A 


123.541 










T-TC33909 


Chrl6 


D 


82.967 










T-TC41731 


Chr20 


D 


80.438 


miR172-Me49b* 


ChrlO 


A 


1 34.949 


T-TC30742 


ChrlO 


A 


123.541 










T-TC33909 


Chrl6 


D 


82.967 










T-TC41731 


Chr20 


D 


80.438 


nniR172-Em20* 


Chr21 


D 


229.243 


T-TC30742 


ChrlO 


A 


123.541 










T-TC33909 


Chrl6 


D 


82.967 










T-TC41731 


Chr20 


D 


80.438 


nniR393-Me4 


ChrOS 


A 


67.257 


T-TC38219 


Chr08 


A 


65.478 


nniR393-Me42 


Chrl2 


A 


154.801 


T-TC38219 


Chr08 


A 


65.478 


miR393-Ennl 


Chrl6 


D 


44.796 


T-TC38219 


Chr08 


A 


65.478 


nniR395-Me4* 


Chrl9 


D 


1 8.742 


T-CO 128820 


Chrl8 


D 


65.452 










T-TC30287 


Chr26 


D 


162.675 


miR396-Me7* 


Chr07 


A 


85.357 


T-TC35055bA 


Chr05 


A 


155.907 










T-TC39199A 


Chrl2 


A 


35.006 










T-TC35055aA 


Chrl9 


D 


95.301 


nniR396-Me43* 


ChrlO 


A 


53.687 


T-TC35055bA 


Chr05 


A 


155.907 










T-TC39199A 


Chrl2 


A 


35.006 










T-TC35055aA 


Chrl9 


D 


95.301 


nniR396a-DX562190* 


Chrl6 


D 


73.216 


T-TC35055bA 


Chr05 


A 


155.907 










T-TC39199A 


Chrl2 


A 


35.006 










T-TC35055aA 


Chrl9 


D 


95.301 


miR397-Me49* 


Chr09 


A 


145.396 


T-AW730829 


Chrll 


A 


226.873 










T-TC31021 A 


Chrl9 


D 


0 


miR397a-DW490000* 


Chr04 


A 


59.354 


T-TC31021 A 


Chrl9 


D 


0 


miR403-Me35 


Chr02 


A 


71.74 


T-CO078480 


Chrl3 


A 


194.48 


miR403-Em27 


Chrll 


A 


146.949 


T-CO078480 


Chrl3 


A 


194.48 


nniR403-Me20 


Chrl7 


D 


84.782 


T-CO078480 


Chrl3 


A 


194.48 


miR479-ES809290a 


Chr26 


D 


44.775 


T-TC31154 


Chrl8 


D 


121.25 



*: miRNA was associated with more than one target on the interspecific BCi linkage map. 

▲ : Target was regulated by more than one miRNAs family on the interspecific BCi linkage map. 



families (Additional file 7: Figure S2L) have the same 
targets as depicted in the network diagram: TC35055 and 
TC39199. Chromosomes on which miR396 and miR396a 
are located are also on homologous chromosomes (Chr07 
and Chrl6), and their target, TC35055, has two sites. Fur- 
thermore, TC350055a, located on Chrl9, and TC350055b, 
located on ChrOS, are found on homologous chro- 
mosomes. Additional file 7: Figure S2G depicts the same 
information: miR196b-Me45 was located on ChrOS and 
miR196b-EmlO was located on Chrl9. For miR156c- 
DX53560 (Additional file 7: Figure S2A), its target, 
TC37339 was located on Chr21 and another target 
TC29212 was located on Chrll, which are homologous 
chromosomes. Additional file 7: Figure S2F depicts the 



location of miR169a (Chrl8) and TC29763 (Chrl3), 
which are also located on homologous chromosomes. 
Other miRNAs and their target genes were located on 
different chromosomes. 

In our network diagram, the relationship and distribu- 
tion of miRNAs and their putative target genes on the 
interspecific BCi genetic linkage map are clearly and in- 
tuitively revealed. This diagram will facilitate our under- 
standing of the regulatory mechanisms and networks for 
target genes and cell development in cotton. Allen et al. 
[36] stated that plant miRNAs are thought to be derived 
from their target sequences after gene duplication, in- 
verted duplication, and divergence in miRNA evolution. 
Additional file 7: Figure S2 shows that many miRNAs and 
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Figure 3 RT-PCR analysis of pre-miRNAs between mapping parents. Numbers at the top represent ODPA, 5DPA, lODPA, 15DPA, 20DPA and 
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Emian22 and 3-79 were classified accordingly. Gene primers are labeled on the left. 







their targets, miRNAs with the same target, and targets 
of the same miRNA family were located on homologous 
chromosomes or the same chromosome. These data are 
consistent with Aliens published findings [36]. Other 
miRNAs and their targets did not reside on the same or 
homologous chromosomes, which may be due to evolu- 
tionary mutations or deletions or insufficient marker num- 
bers. It was noteworthy that homologous chromosomes, 
Chr5 and Chrl9, contained more targets and mirRNAs as 
depicted in the network diagram, and this may have impli- 
cations for future studies of miRNAs and their targets. 

Expression differences between parents by RT-PCR 
and qRT-PCR 

We used RT-PCR analysis to investigate pre-miRNAs 
and their putative target genes during fiber development 
between G. hirsutum and G. harhadense, A total of 181 
primer pairs were used for RT-PCR analysis, including 
83 pre-miRNA primers and 98 randomly selected target 
gene primers. Twenty- five pre-miRNA genes (30.1%) 
were expressed in cotton fibers. Among them, 88% of 
these expressed pre-miRNAs were significantly different 
among various stages of fiber development (0, 5, 10, 15, 
20, and 25 DPA) between Emian22 and 3-79. Ten pre- 
miRNAs had weak or no expression in both Emian22 
and 3-79 (defined as no expression differences in this 
study). Furthermore, 13 pre-miRNAs had similar ex- 
pression patterns, and 9 pre-miRNAs had different 



expression patterns (Figure 3). Figure 3 also shows that 
40% of the pre-miRNA genes had higher expression in 
early ovary development. 

Regarding target genes, 66 (67.3%) were expressed in 
cotton fiber (Additional file 8: Table S6). Among them, 50 
were obviously differentially expressed; 14 had minor dif- 
ferences; and 2 were not different (Figure 4). Overall, 
75.8% of the expressed target genes were significantly dif- 
ferent at various fiber development stages (0, 5, 10, 15, 20, 
and 25 DPA), when comparing Emian22 and 3-79. Three 
target genes were weakly expressed or not expressed at all 
stages in both Emian22 and 3-79 (defined as no expres- 
sion differences in this study). Furthermore, 39 target 
genes had similar expression patterns. Of the 25 target 
genes with different expression patterns, 13 were up- 
regulated and 12 were down regulated in 3-79, compared 
with Emian22. To further confirm RT-PCR data, target 
genes belonging to different categories were randomly 
chosen for qRT-PCR analysis (Additional file 9: Figure 
S3). Consistent results were observed between both RT- 
PCR and qRT-PCR analyses. 

In this study, a diagram was constructed to depict 
pre-miRNAs and their corresponding target genes 
(Additional file 10: Figure S4). Additional file 10: Figure 
S4 shows that miR156c-DX535560, miR472-DT527030, 
and miR482a-DR457519 were similar with respect to 
their targets in various fiber development stages of 
Emian22. However, no correlation was observed between 
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Figure 4 RT-PCR analysis of target genes between mapping parents. Numbers at the top represent ODPA, 5DPA, lODPA, 15DPA, 20DPA and 
25DPA of Ennian22 and 3-79, respectively. Similar expression tendencies between Emian22 and 3-79 were classified into similar expression 
patterns. Obvious differential expression tendencies between Emian22 and 3-79 were classified into differential expression patterns. Obvious 
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Emian22 and 3-79 were classified accordingly. Gene primers and their corresponding miRNA families are labeled on the left, and only one miRNA 
family is listed. More family information was shown in Additional file 8. 



pre-miRNA and its target expression in various fiber de- 
velopment stages of 3-79. MiR157-DT558905, miR395a- 
DT567568 and miR395a-DT567568 had contrasting 
expression tendencies with their targets in various sta- 
ges of Emian22, but no correlation was observed be- 
tween the expression of pre-miRNA and its targets 
in various fiber development stages of 3-79. MiR172- 
TC124126, miR413-AW187128, and miR414a-DR452834 
and their targets had similar expression tendencies in 



various parent stages. MiR414d-DW514431 had contrast- 
ing expression tendencies with its targets in various stages 
of Emian22 and had similar expression tendencies in vari- 
ous stages of 3-79. 

In summary, we report dynamic expression of pre- 
miRNAs and putative target genes at various fiber develop- 
ment stages in G. hirsutum and G. barbadense. Most 
expressed pre-miRNAs (88%) and their putative targets 
(75.8%) were obviously different with respect to expressions 
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at different stages between the two parents. Unique expres- 
sion patterns of pre-miRNAs and their putative targets may 
be connected with a particular function. We also found that 
pre-miRNA families and their putative targets had correlated 
expression tendencies in various growth stages, suggesting 
that miRNA genes may regulate target gene expression. 
However, the biological relevance of this must be investigated 
further to better understand regulatory mechanisms and the 
overall network of plant growth control. 

Conclusions 

In this study, miRNAs and target markers were developed. 
We found that although the primers were specific, pre- 
miRNA primer polymorphisms were low (14.5%), and tar- 
get primer polymorphisms were even lower (11.7%). The 
63 miRNA loci were evenly distributed on 24 chromo- 
somes, except Chr22 and Chr24. Also, 156 target gene loci 
were preferentially distributed on the Dt sub-genome and 
some chromosomes. Genetic mapping in a segregating 
population provides an opportunity to examine distribu- 
tion of miRNAs and their target genes. A network dia- 
gram depicting miRNAs and their targets confirmed that 
plant miRNAs may be derived from their target sequences 
after gene duplication or inverted duplication [36]. Com- 
parative expression analyses of Gossypium hirsutum and 
G. barbadense revealed that miRNAs and their target 
genes play a role in cotton fiber development. 

Methods 

Marker development 

Pre-miRNA sequences reported in the literature and a 
database (http://www.mirbase.org) were downloaded 
[37-43] (Additional file 1: Table SI). Putative target genes 
were predicted and downloaded with psRNATarget (http:// 
plantgrn.noble.org/psRNATarget/). Because one target gene 
may relate to more than one miRNA, repeat target se- 
quences were deleted. Finally, a total of 123 pre-miRNA 
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sequences and 1,399 putative target sequences were ob- 
tained. Primers were then designed based on the pre- 
miRNAs and the target sites of the putative target genes 
by Primer 3.0 (examples in Figure 5). Criteria for primer 
design were a primer length of 18-25 bp (20 bp 
optimum), GC content of 35-70% (50% optimum), 
annealing temperature of 50-65°C (55°C optimum), and 
PGR product size of 100-1,000 bp. 

Because the polymorphism is low in cotton, to genetic- 
ally map more miRNAs, SRAP primers, including 64 
forward primers and 64 reverse primers [44] were com- 
bined with 23 miRNA degenerate primers [32] in the 
form of miRNA-SRAP analysis. Finally, a total of 2,944 
primer combinations were applied in this study. 

PGR amplification and electrophoresis 

PGR of pre-miRNAs and their putative target genes was 
performed in a solution (10 (iL) containing 25 ng DNA 
template, 1 x Buffer, 2.0 mmol L'^ MgGl2, 0.25 mmol L'^ 
dNTPs, 0.16 (imol L'^ of forward primer, 0.16 (imol L'^ 
of reverse primer, 0.8 units of Taq DNA polymerase, and 
ddH20 was added to 10 (iL. The PGR program was 
performed using the following profile: 95°G for 5 min, 
followed by 34 cycles consisting of 94°G for 50 sec, 56°G 
for 45 sec, and 72°G for 60 sec; and a final extension 
step of 5 min at 72°G. Then, PGR products were sepa- 
rated on an 8% non-denaturing polyacrylamide gel at a 
constant voltage (15 W) for about 4 h at room 
temperature. After electrophoresis, DNA fragments were 
detected by silver staining, coloring in a sodium hydrox- 
ide and formaldehyde solution. 

PGR of an miRNA-based SRAP marker was per- 
formed in a solution (10 (iL) containing 30 ng DNA 
template, 1 x Buffer, 2.0 mmol L'^ MgG^, 0.2 mmol L'^ 
dNTPs, 0.1 (imol L'^ of forward primer, 0.1 (imol L'^ of re- 
verse primer, 0.5 units of Taq DNA polymerase, and ddH20 
was added to 10 (iL. The PGR program was performed 
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Figure 5 Primer design strategies for miRNAs and target genes, (a) miR414d from DW5 14431 (Genbank acc No. DW5 14431); (b) target gene 
of miR414d from TC31963 (Cotton Gene Index No. TC31963, http://compbio.dfci.harvard.edu/tgi/cgi-bin/tgi/gimain.pl?gudb=cotton). *: Target 
sequence; >>>>>: Primer region. 
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using the following profile: 94°C for 5 min, followed by 4 
cycles of 94°C for 1 min, 35°C for 1 min, and 72°C for 1 
min; then by 34 cycles of 94°C for 1 min, 50°C for 1 min, 
and 72°C for 1 min; and a final extension step of 10 
min at 72°C. Then, PGR products were separated on a 
6% denaturing polyacrylamide gel (29:1 acrylamide and N, 
N-methylene bisacrylamide) at a constant voltage (80 W) 
for about 2 h at room temperature. After electrophoresis, 
DNA fragments were detected by silver staining, coloring 
in a sodium hydroxide and formaldehyde solution. 

Genetic mapping 

Our BCi mapping population ([Emian22 x 3-79] x 
Emian22) [45,46] was used to map miRNAs and their 
putative target genes. All primers with polymorphisms 
between parents were used to enrich the interspecific 
BCi genetic linkage map [46]. Polymorphic loci were 
integrated into the interspecific BCi linkage map using 
JoinMap V3.0 [47]. The logarithm of odds (LOD) 
threshold was 5.0. Map distances in centi Morgans 
(cM) were calculated using the Kosambi mapping func- 
tion [48]. The linkage map was drawn by MapChart 
V2.2 software [49]. 

Drawing method of network diagram 

In this study, miRNAs and their putative target genes were 
mapped in cotton chromosomes. The distribution of 
miRNAs and their putative target markers was obtained 
from the genetic linkage map. MiRNA loci and their puta- 
tive target loci that were mapped on the chromosomes 
were then "connected" using individually colored lines to 
indicate each miRNA family's connection to their putative 
target genes. Then, a network diagram was obtained based 
on the distribution of miRNAs and their target markers 
on genetic linkage map by Circos (http://circos.ca/). 

RT-PCR and qRT-PCR analysis 

To determine the expression difference of pre-miRNAs 
and their putative target genes during fiber development 
between G. hirsutum and G. harhadense, RNAs were 
extracted from developing fibers at 0 days post anthesis 
(DPA), 5 DPA, 10 DPA, 15 DPA, 20 DPA and 25 DPA, 
and RNAs (4 (ig) were reverse-transcribed to cDNA by 
the M-MLV-RT Reverse Transcriptase (Invitrogen). RT- 
PCR and qRT-PCR analyses were performed according to 
the methods described by Munis et al. [50] with minor 
modifications. Ubiquitin (GenBank acc No.: DQl 16441) 
primer pair (forward 5'GAAGGCATTCCACCTGACCA 
AC3', reverse 5 'CTTGACCTTCTTCTTCTTGTGCTTG 
3') [51] was used as an internal standard showing the 
equal amounts of first-stranded cDNA in each sample. In 
addition, qRT-PCR experiments were biologically repeated 
three times. 



Additional files 



Additional file 1: Table SI. Source of miRNA sequences and 
distribution in the scientific literature. 

Additional file 2: Table S2. Characteristics of 83 pre-miRNA primers in 
this study. V: Polymorphic primers; x: No polymorphic primers: The 
primer was not mapped on the chromosome. 

Additional file 3: Table S3. Characteristics of 1,255 target primers in 
this study. V: Polymorphic primers; x: No polymorphic primers. -: The 
primer was not mapped on the chromosome. 

Additional file 4: Table S4. Characteristics of 55 miRNA-SRAP 
polymorphic primer combinations in this study-: The primer was not 
mapped on the chromosome. 

Additional file 5: Figure SI. Locations of polymorphic markers on the 
BCi genetic linkage map. miRNA markers are underlined and in bold. Target 
markers are underlined. miRNA-based SRAP markers are italic and in bold. 

Additional file 6: Table S5. Marker distribution on the interspecific BCi 
genetic linkage map. 

Additional file 7: Figure S2. Network diagram revealed by 
relationships between one miRNA family and their target genes. Families 
are as follows: A: miR156c family; B: miR157b family; C: miR159 family; D: 
miR164 family; E: miR167 family; F: miR169a family; G: miR169b family; H: 
miR172 family; I: miR393 family; J: miR395 family; K: miR396 family; L: 
miR396a family; M: miR397 family; N: miR397a family; 0: miR403 family; P: 
miR479 family. Red circle: MiRNAs; black circle: Target genes. The scale 
marked on each chromosome represented a genetic map distance (cM). 

Additional file 8: Table S6. Target gene primers and their 
corresponding miRNA families. 

Additional file 9: Figure S3. qRT-PCR analysis of eight primers 
randomly chosen from pre-miRNAs and targets. Expression levels of 
Emian22 and 3-79 are shown. represents P<0.05, and represents 
P<0.01. Primers are as follows: A: miR164-BF270853; B: miR393-TCl 18931; 
C: miR413-AW187128; D: T-BF275332; E: T-BM3579942; 
F: T-TC27422; G: T-TC37575; H: T-TC38767. 

Additional file 10: Figure S4. Wiring diagram of target gene primers 
and their corresponding miRNA families, (a) RT-PCR analysis of miRNA 
genes; (b) RT-PCR analysis of target genes. 
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